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We study the effects of next-nearest neighbor (NNN) interactions in the two-dimensional fer- 
romagnetic kinetic Ising model exposed to an oscillating field. By tuning the interaction ratio 
(p — Jnnn / Jnn) of the NNN (Jnnn) to the nearest-neighbor (NN) interaction (Jnn) we find that 
the model undergoes a transition from a regime in which the dynamic order parameter Q is equal 
to zero to a phase in which Q is not equal to zero. From our studies we conclude that the model can 
exhibit an interaction induced transition from a deterministic to a stochastic state. Furthermore, we 
demonstrate that the systems' met ast able lifetime is sensitive not only to the lattice size, external 
field amplitude, and temperature (as found in earlier studies) but also to additional interactions 
present in the system. 
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I. INTRODUCTION 

The two-dimensional (2D) ferromagnetic nearest- 
neighbor kinetic Ising (NNKI) model in an oscillating 
field has been used extensively to study non-equilibrium 
(NEQ) properties [BE]- The presence of an oscillatory 
field introduces an explicit time dependence in the Hamil- 
tonian. This causes the system to exhibit a hysteretic re- 
sponse. Below the critical temperature (Tc) the system 
exhibits a dynamic phase transition (DPT) or stochas- 
tic resonance (SR) [l]-[5] based upon the strength of the 
field amplitude (ho)^ frequency (/), temperature (T) of 
the system, and the lattice size (L). These parameters 
dictate whether the system will be in a deterministic or a 
stochastic regime. These regimes are further subdivided 
into one of the four possible regions: strong field (SF) 
and multi-droplet (MD) for deterministic, single-droplet 
(SD) and co-existence (CE) for stochastic. The SF and 
the MD regions exist for strong fields/large system sizes 
and the SD exists for weak fields/small system sizes [2]. 
The MD and SD regions have been studied exhaustively 
for metastability mechanisms, finite-size scaling effects, 
DPT, SR, hysteresis exponents, universality, critical ex- 
ponents, and effect of square-wave oscillating external 
field with a soft Glauber dynamics [3H8l . 

A measure of the DPT is given by the dynamic or- 
der parameter, Q, which is the period- averaged magne- 
tization. The DPT occurs between an ordered dynamic 
phase with {\Q\) ^ and a disordered dynamic phase 
with (IQI) = [2 (only in the MD region). The DPT 
in the NNKI model is of second order [2]-[4]. Recently, 
evidence for a dynamic phase transition has been inves- 
tigated experimentally in [Co/Ptja magnetic multilayers 
|9 . Furthermore, Monte Carlo simulations and mean- 
field studies of hysteresis in the NNKI model show that 
the loop area undergoes a transition from a symmetric 
shape to an asymmetric shape or vice- versa [1, 10 . The 
breakdown in the shape of the hysteresis loop has also 
been attributed to a DPT with a spontaneously broken 



symmetric phase [T]. The scaling relations for the hys- 
teresis loop area have been measured in ultrathin and 
thin ferromagnetic film systems such as Fe/Au(001) [11], 
FesoNigo [12 , and Co/Cu(001) [13 . In the SR region 
the magnetization switches through random nulcleation 
of a single droplet of spins aligned with the applied field. 
This region has been studied for stochastic hysteresis us- 
ing time-dependent nucleation theory and a variable rate 
Markov processes [6 . 

In this paper we extend the NNKI model to include 
next-nearest neighbor (NNN) feromagnetic interactions. 
We term this model the next-nearest neighbor kinetic 
Ising (NNNKI) model (see Eq. [T]). The motivation be- 
hind this investigation is to study the effects of addi- 
tional (NNN) interactions in the system. In Sec. [ll] we 
state the NNNKI model and describe the Monte Carlo 
method used. In Sec. |III| we present our results on the 
NNNKI model. In SecTTlVl we discuss the results and 
state the main conclusions of our paper. 

II. MODEL AND METHOD 

The model used in this study is the kinetic NNN Ising 
ferromagnet on a square lattice with periodic boundary 
conditions. The NNNKI Hamiltonian is given by 

H = - Jnn ^ SiSj - Jnnn ^ SiSj - ho sin(27r/t) ^ 5^,(1) 

{id) [id] i 

where Si is the ith spin and can have values of Si = ± 
1, Jnn is the NN coupling, Jnnn is the NNN coupling, 
ho is the external field amplitude, and / is the frequency 
of the external field. The sums Xl(^j) and "^[i j] run 
over all NN and NNN pairs, respectively. Both the cou- 
plings are ferromagnetic, Jnn > and Jnnn > 0- The 
ratio of the couplings is defined as p = Jnnn/ Jnn- The 
spin-flip dynamics used is the Metropolis algorithm with 
the Monte Carlo step per spin (MCSS) as the unit time 
step [T4^. The system is allowed to be in contact with a 
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heat bath at temperature T, and each attempted spin 
flip from Si —Si is accepted with the probability 
W{Si -Si) = exp(-/3A E^). Here A^;^ is the change 
in energy of the system that would result if the spin flip 
were accepted and (3 = l/Zc^T where the Boltzmann con- 
stant /cb is set equal to one. Using the above Hamilto- 
nian (Eq. [T]) and the Monte Carlo method we compute 
the dynamic order parameter Q 

Q = ^W^^^, (2) 

I 



where m{t) is the time dependent magnetization per unit 
site. In the next two sections we study the physics of the 
NNNKI model. 




FIG. 1: For L = 128, K = 0.5Jnn (where Jnn = 1), / = 10"^ 
and T — 0.8T^^ we have (a) Magnetization time series data 
demonstrating the deterministic to stochastic transition for 
interaction ratio p = - 0.4. The solid line represents the 
systems magnetization and the dashed line is the external 
field, (b) Spontaneous symmetry breaking of the hysteresis 
loops for the same set of interaction ratios as in part (a). 



III. RESULTS 

A. Interactions and their effects 

In all the previous studies involving the NNKI model 
the frequency, magnetic field or temperature was used 
as the tuning parameter to demonstrate a transition be- 
tween the stochastic to deterministic region. This al- 
lowed the study of SR or DPT based upon the system 
size and field amplitude. Our goal in this section is to 
demonstrate that the deterministic to stochastic tran- 
sition can also be achieved by tuning the interactions. 
In our computation we choose a set of (ho^f^T) values 
so that the system is in the deterministic MD region in 
the NN model. We take ho = 0.5Jnn, / = 10"^, and 
T = 0.8T^^. A choice of (/io,f,T) so that the system is 
in stochastic region is also equally valid. We then change 
the values of interaction ratio, p, and compute the mag- 
netization time series data for several values of p (see 
Fig. [1(a)] ). In our calculations we take p to range from 
to 0.4. This range of parameters is sufficient to high- 



light the physics of the problem. From the magnetization 
time series data we see that the system makes a transition 
from a deterministic to a stochastic region. For the same 
values of p we plot the hysteresis loops which show the 
spontaneous symmetry breaking (see Fig. |l(b)" ). We also 
track the values of the dynamic order parameter {\Q\) 
(see Fig. [2| and we find that the value of {\Q\) changes 
from zero to a non-zero number with increasing p. The 
system is therefore sensitive to NNN interactions. This 
sensitivity can be physically explained by considering the 
effect of interactions on the metastable lifetime. 



B. Metastable lifetime, r{ho,T, L,p) 

The key to understanding the observed phenomena in 
the NNNKI model is to study the metastable lifetime. To 
determine r for this model we performed several instan- 
taneous field reveral simulations p^. The system was 
initially prepared in an all up spin configuration. The 
field was instantaneously reversed and the relaxation of 
the system then studied. The metastable lifetime is the 
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FIG. 2: Tracking the dynamic order parameter, (|Q|), as the interaction ratio, p = Jnnn/Jnn, is changed. The value of (|Q|) 
changes from a zero to a non-zero number. The lattice sizes are L=64, 72, 96, and 128. The computations were done at 
ho = O.bJnn (where Jnn = 1), / = 10~^, and T = 0.8T^^ . Each run is of length 100,000 MCSS. The dashed lines are a guide 
to the eye. 



number of MCSS needed for the system to decay to a 
net zero magnetized state {m{t) = 0) from a completely 
magnetized state (m(t) = ±1). The average metastable 
lifetime, {r{ho^T^ L^p)) ^ is calculated after 1000 repeated 
trials. The simulations were performed for L=128 at 
T = O.ST^^ for various ratios of the interaction strength 
p=0, 0.5, 0.7, and 1. The results are shown in Fig. |3] 
From Fig. [3] we see that the NNNKI model's average 
(r(/io, T, is much greater than the NNKI model's 

(r(/io, T, The lifetime is sensitive to the ratio of 

the interaction strengths. Since the metastable lifetime 
dictates the underlying metastable mode decay this in 
turn causes the system to exhibit a transition from the 
deterministic to the stochastic region. 
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IV. DISCUSSION AND CONCLUSION 



The computations in this article demonstrate that in- 
teractions have an important role to play in the kinetic 
Ising model in an oscillating applied field. While we have 
established that the system undergoes a transition from 
a deterministic to a stochastic regime our computations 
do not highlight the DPT or SR possible in the system. 
Some of the magnetization time series data (p = 0.22 
to 0.25) show an inkling of a MD region as seen in the 
"wandering" time series graphs. However, the response 
of the system eventually becomes stochastic. Our choice 
of a small lattice size (L = 128) and system parame- 
ters (I/,/io, and T) may have caused the MD region to 
be less prominent compared to SR since the SR phe- 
nomenon is a finite-size effect. Investigating the model 
further for a DPT and/or SR for changing interaction 
strengths will be a topic of future investigation to be re- 
ported elsewhere. Another important conclusion of this 
article is that the metastable lifetime, r, which dictates 
the physics of the kinetic Ising model is sensitive to addi- 
tional interactions. Finally, it is our hope that this work 
will motivate experimentalists to investigate the interac- 
tion induced transition observed in the NNNKI model in 
real material systems. 



FIG. 3: Variation of the average metastable lifetime, 
{T(ho,T, L,p)) , as a function of interaction strength ratio 
P = Jnnn/Jnn ou a log-llucar scalc. The average lifetime 
was computed using the first passage of time to zero mag- 
netization for 1000 repeated trials for different values of the 
external magnetic field amplitude ho. The calculations were 
performed for L=128 at T = 0.8T^^ . We see from the plot 
that the value of (r) is sensitive to additional interactions 
present in the system. The computed error bars (standard de- 
viation) are displayed. For some of the computed data points 
the error bars are smaller than the symbol size. 
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